Highly charged ions with El, Ml, and E2 transitions within laser range 
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Level crossings in the ground state of ions occur when the nuclear charge Z and ion charge Zi on 
are varied along an isoelectronic sequence until the two outermost shells are nearly degenerate. 
We examine all available level crossings in the periodic table for both near neutral ions and highly 
charged ions (HCIs) . Normal El transitions in HCIs are in X-ray range, however level crossings allow 
for optical electromagnetic transitions that could form the reference transition for high accuracy 
atomic clocks. Optical El (due to configuration mixing), Ml and E2 transitions are available in 
HCIs near level crossings. We present scaling laws for energies and amplitudes that allow us to make 
simple estimates of systematic effects of relevance to atomic clocks. HCI clocks could have some 
advantages over existing optical clocks because certain systematic effects are reduced, for example 
they can have much smaller thermal shifts. Other effects such as fine-structure and hyperfine 
splitting are much larger in HCIs, which can allow for richer spectra. HCIs are excellent candidates 
for probing variations in the fine-structure constant, a, in atomic systems as there are transitions 
with the highest sensitivity to a- variation. 
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I. INTRODUCTION 

Current technological plans hint to mainstream adop- 
tion of highly charged ions (HCIs) for many uses in the 
near future (see, e.g. the review [l|). Production of any 
ion stage of practically any naturally occurring element 
is possible at ion accelerators and/or electron beam ion 
traps. Furthermore great progress has been made re- 
cently in trapping, cooling, and spectroscopy of HCIs 
(see, e.g. and the review @). In this paper we con- 
sider candidate transitions for an optical clock made us- 
ing a HCI that has a configuration crossing in the ground 
state: a "level crossing" . 

Level crossings in ions occur when the energy order- 
ing of orbitals changes with increasing ion charge. The 
ion charge may be increased by considering ionisation 
along an isonuclear sequence or by considering an isoelec- 
tronic sequence with variable nuclear charge. The latter 
is somewhat simpler to deal with theoretically since the 
electronic structure does not usually change very much 
between adjacent ions. In this paper we discuss isoelec- 
tronic sequences at points where the electronic struc- 
ture does change — the level crossings — and interesting 
properties can emerge. Near level crossings, the frequen- 
cies of transitions involving the crossing orbitals can be 
much smaller than the ionisation energy. This means 
that they can be within the optical range and have the 
potential to be excited by lasers, opening the possibility 
of performing high-precision spectroscopy and building 
optical clocks using HCI reference transitions. 

This work is also motivated by astronomical observa- 
tions of quasar absorption spectra that hint that there 
is a spatial gradient in the value of the fine-structure 
constant, a = e 2 /he 0, S|. Data samples from the 
Very Large Telescope and Keck Telescope 0, H3] inde- 
pendently agree on the direction and the magnitude of 
this gradient, which is significant at a 4.2<r level. A con- 
sequence of the astronomical result is that since the solar 
system is moving along this spatial gradient, there may 



exist a corresponding temporal shift in a in the Earth 
frame at the level a/a ~ 10~ 19 yr _1 [ll|. Finding this 
variation using atomic clocks could independently cor- 
roborate the astronomical result in the laboratory. 

The best current terrestrial limit on time- variation of 
a was obtained by comparing the ratio of frequencies of 
the Al + clock and the Hg + clock over the course of a 
year [l2j]. The ratio is sensitive to a- variation because 
the reference transitions in the two clocks have different 
sensitivity coefficients, q, defined as 
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dx 
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where x — o? /a\ — 1 is a normalised change in a 2 from 
the current value a\ , and q and tu are measured in atomic 
units of energy. In this experiment the Al + clock is 
relatively insensitive to a variation (low q coefficient), 
thus serving as an "anchor" line. On the other hand the 
Hg + clock is sensitive to a- variation (high q coefficient). 
Therefore, the ratio of these transition frequencies will 
change if a changes. The limit on the rate of change of 
a was measured as a/a = (—1.6 ± 2.3) x 10~ 17 yr _1 . 

To compete with astrophysical measurements of the 
spatial gradient, the atomic clock limits must be im- 
proved by around two orders of magnitude. Several pro- 
posals have been made for atomic clocks that, if measured 
at the same level of accuracy as the Al + /Hg + ratio, would 
give much stronger limits on a- variation. These include: 
proposals to construct clocks using heavier elements with 
similar properties (e.g. the Tl + clock proposed by [HI); 
systems with large relative sensitivities to a- variation ex- 
ploiting the accidentally degenerate levels in Dv [Til HH 
or fine-structure anomalies in Te, Po, and Ce [16j ; a va- 
riety of transitions in heavy elements with large q val- 
ues (e.g. [T7I - I20I ): and nuclear clocks based on the 7.6eV 
isomeric transition in the 229 Th nucleus that would have 
extraordin ary sensitivity to variation of fundamental con- 
stants (21H24 | . For a more complete review see [2^, [26[ . 

Transitions near level crossings in HCIs can provide 
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higher sensitivity to a-variation than any other optical 
transitions seen in atomic systems [27], [H[. Consider the 
following analytical formula for the relativistic shift of an 
energy level in the single-particle approximation [l4j |: 



(Zaf 



(2) 



where /„ is the ionisation energy of the orbital (atomic 
units h = e = m e = 1), and v is the effective principal 
quantum number. A transition in a HCI can have a large 
sensitivity because the difference in q n between the lev- 
els involved can be large. The enhancement comes from 
the coherent contributions of three factors: high nuclear 
charge Z, high ionisation degree Z loll (leading to large 
7 n ), and significant differences in the configuration com- 
position of the states involved (large changes in j and 
v). For nearly- filled shells, an additional enhancement in 
the a-sensitivity occurs due to each electron spending ap- 
proximately half its time nearer to the nucleus than other 
electrons in the same shell. In these cases q n ~ ifj 2 [28[ . 

In this paper we perform a systematic search for level 
crossings in HCIs throughout the periodic table. We 
identify several ranges of Z and Z ion where level cross- 
ings can be found, and perform configuration interaction 
calculations for some of the most promising systems. In 
Sec. lVIII we discuss how systematic effects that affect op- 
tical clocks arc modified in the case of HCIs, and find that 
HCIs confer some benefits over near-neutral ions. Cur- 
rent experimental techniques might be applied to build 
a similar clock retaining high precision, but with much 
higher sensitivity to a- variation. 



II. METHOD 

Our first task in this work is to identify HCIs with level 
crossings in the ground state. We start with neutral ions 
and then increase Z, working along the isoelectronic se- 
quence from the neutral atom filling order towards the 
Coulomb filling order. The Madelung rule (also known 
as the Klechkowski rule) can be taken as a first approxi- 
mation for determining the filling order of electron shells 
in neutral atoms. We show in the Appendix that this is a 
good approximation because deviations from this filling 
order in neutral atoms disappear with a small increase in 
the ion charge Z- lcm . Also we know that in very highly- 
charged ions, the energy levels of the electrons must ap- 
proach the hydrogen-like (Coulomb) limit, where all or- 
bitals with the same principal quantum number n are 
nearly degenerate. Fig. [TJ presents the order of electron 
orbitals under both ordering schemes. Since any differ- 
ence in the ordering as computed from the Madelung 
rule and that of the hydrogen-like limit must be resolved 
with increasing Z ion , the 'out-of-order' levels must cross 
at some Z lon . From the transition between these limits it 
is seen that the only types of crossings available in HCIs 
are between orbitals with angular momenta s — d, s — /, 
and p — f. 



FIG. 1. A comparison of the ordering of electron orbitals: 
The first column is the order of filling as derived by applying 
the Madelung rule, while the second column is derived for 
a hydrogen-like atom (excluding gr-wave and /i-wave orbitals 
that cannot be occupied in the ground state of any real ion) . 
The ordering of orbitals changes with increasing ion charge, 

Zion ■ 

Madelung Ordering Coulomb Ordering 



Is 

2s 

2p 

3 s 

3p 
4s 
3d 
4p 
5 s 
4d 
5p 
6s 
4f 
5d 
6p 
7s 
5f 
6d 




Is 

2s 

2p 

3s 

3p 
3d 
4s 

4p 
4d 
4f 

5s 

5p 
5d 
5f 
6s 
6p 
6d 
6f 
7s 



Neutral atoms sometimes have ground state electronic 
configurations that deviate from the Madelung rule. In 
isoelectronic sequences starting with such atoms other 
types of level crossings can occur (namely, 5d — 4/ and 
6c? — 5/). However, we find that the new crossings occur 
with the addition of just a few extra protons; no addi- 
tional crossings are found in highly- charged (Z lon > 5) 
ions. Full details are presented in the Appendix. 

To find the ions along the isoelectronic sequence where 
level crossings lead to small transition frequencies we 
perform Dirac-Fock (relativistic Hartree-Fock) calcula- 
tions. An example is presented in Fig. [5J which shows 
the 4/ and hp valence orbitals of the indium isoelec- 
tronic sequence (N = 49) calculated using Dirac-Fock 
in the V 1 ^^ 1 approximation. It is seen that for low val- 
ues of Z- lon the 4/ orbitals lie above the hp orbitals, but 
at Z = 59 the 4/ levels drop below the hp 3 / 2 orbital, 
and between Z — 59 and 60 they cross the 5pi/2 orbital 
energy. In general, this method produces acceptable es- 
timates for the position of the crossings as we will see by 
comparison with Configuration Interaction calculations 
in Sections fVl and fVTl 

Many-body perturbation theory (MBPT) corrections 
can be included, but our calculations show that this does 
not change the position of the crossing point (see Fig. [31 
which shows a detailed view of the crossing point in Fig. [5] 
with and without including of MBPT corrections). 
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FIG. 2. Dirac-Fock energies of the 4/ 5 / 2 (solid), 5pi/2 
(dashed), and 5p3/2 (dotted) levels of the In (N = 49) iso- 
electronic sequence. 
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The indium sequence described above has one valence 
electron above closed shells (a cadmium core that we may 
consider frozen). In general, however, we can perform 
Dirac-Fock calculations even for several- valence-electron 
ions provided we scale the contribution of each subshell 
by its filling fraction. Again this gives reasonable ac- 
curacy for the ionisation energy (order of few percent) 
which is good enough to identify level crossings. As we 
progress along an isoelectronic sequence, we increase Z 
until the first crossing point is reached. After this point 
the electronic configuration will be altered, and to find 
other crossing points that occur later in the sequence 
further calculations must be performed with the modi- 
fied electron configuration which assumes that the first 
crossing has occurred. 

In principle, it is possible to use the weighted Dirac 
Fock method outlined above for an arbitrary number of 
electrons, but for partially-filled shells and electron-hole 
calculations there usually will be more than one possible 
DF electron configuration to use. One such example is 
Cr II [2^] , where the d-shell electrons must be accounted 
for in the DF approximation, but it is not clear if a V N 
scheme where 3d 5 is included in the DF potential or a 
V N ~ X scheme where 3d 4 is included will give better agree- 
ment with experiment (of course, in the limit of a com- 
plete basis set both approximations will give the same 
CI result). Furthermore, using a poor approximation for 
the DF potential may result in the Dirac-Fock calculation 
showing no available level crossings. In order to resolve 
this issue, we must perform at least minimal configura- 
tion interaction calculations to locate the crossing point 
and calculate approximate transition frequencies. 

HCIs with many valence electrons have some benefits 
for potential clock applications because of the availability 
of different angular momentum states and configurations. 
This is useful both for finding reference transitions with 
desirable properties and also for increasing the sensitivity 
of the transition to cv-variation, as the q values for an k- 



FIG. 3. Energies of the 4/ 5 / 2 (solid), 5pi/ 2 (dashed), and 
5^3/2 (dotted) levels of the In (N = 49) isoelectronic sequence. 
Upper panel: detail of the level crossing in the Dirac-Fock 
approximation of Fig. [2] Lower panel: the same level crossing 
calculated with many-body perturbation theory corrections 
included. The qualitative nature of the crossing point is not 
significantly affected by the MBPT corrections. 
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electron transition is approximately k times the q value 
for the single electron transition. This is illustrated, e.g., 
by the examples presented in [2?], H^- Furthermore, using 
configuration mixing it is also possible to generate El 
transitions using multiple electrons mas — / crossing. In 
the following sections, we will list all the available level 
crossings in elements from the considerations discussed 
above. 

The calculations presented in this paper use the atomic 
structure code AMBiT ^30], which includes Dirac Fock 
(DF) and Configuration Interaction (CI) algorithms. 
While core-valence calculations can be included in a CI 
calculation via many-body perturbation theory using the 
CI+MBPT method [3l|, for our current purposes this 
is not required, as discussed below. All CI calculations 
for two and three- valence-electron ions are performed us- 
ing a fairly small B-spline basis of the type developed 
in (32l. [33|. including valence orbitals only up to 7spdf. 
MBPT corrections are much more important for the cal- 
culation of transition frequencies, us. More precise studies 
of those HCIs that are of interest to experimentalists will 
need to be performed using the full CI+MBPT theory. 

In Table|T]we compare experimental ionization energies 
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TABLE I. Dirac-Fock calculation of ionization energy and en- 
ergy intervals for neutral lithium, compared with available 
experimental data from [33 ]. 



Level 


J 


Ionization Energy ( 


cm 1 ) 


% Deviation 






DF Calc. 


Expt. 




2s 


1/2 


-43087 


-43487 


-0.919 


2v 


1/2 


-28232 


-28583 


-1.226 




3/2 


-28232 


-28583 


-1.227 


3s 


1/2 


-16197 


-16281 


-0.513 


3p 


1/2 


-12460 


-12561 


-0.808 




3/2 


-12459 


-12561 


-0.809 


3d 


3/2 


-12194 


-12204 


-0.079 




5/2 


-12194 


-12204 


-0.079 


4s 


1/2 


-8444 


-8475 


-0.367 


4p 


1/2 


-6975 


-7017 


-0.604 




3/2 


-6974 


-7017 


-0.604 


4d 


3/2 


-6859 


-6863 


-0.065 




5/2 


-6859 


-6863 


-0.065 



with those calculated in the DF approximation for several 
levels of neutral lithium. For this simple case, we see 
that the ionization energies and intervals are accurate 
to ~ 1% or better. Table [TT| compares the ground state 
ionization energies for selected ions along the tungsten 
(Z = 74) isonuclear sequence with available data. In 
HCIs, we see that we maintain roughly the same degree 
of accuracy, therefore, we can surmise that the position 
of level crossings are fairly accurately determined from 
the DF calculations alone. 

On the other hand, we note that the energy of transi- 
tions between the levels participating in the optical level 
crossing may not be as easily determined for HCIs (in- 
deed even the ordering may be difficult to determine). 
This is because we are selecting HCIs where the differ- 
ence between the ionization energies of these levels are 
strongly suppressed. For the typical scale of ionization 
energies in HCIs, ~ 10 7 cm -1 , an interval of ~ 10 4 cm -1 
is the result of a cancelation at the level 99.9%. To deter- 
mine the ground state conclusively an accuracy of better 
than 0.1% in the ionisation energy is required. As a re- 
sult, in HCIs near level crossings the electronic structure 
is not as well determined as in near neutral ions, and the 
ground state is typically not identified to a high degree 
of confidence. For experimental purposes however, the 
two (or more) possible ground states in HCIs with op- 
tical level crossings can all be considered metastable, as 
they have typical lifetimes ranging from seconds to the 
lifetime of the universe. 



III. GROUND STATE LEVEL CROSSINGS 

In this section we list all possible level crossings that 
occur because of the transition from the Madelung filling 
scheme to the Coulomb degenerate scheme. All crossings 
in Fig.Q]are represented, however most occur at relatively 
low ion stage or outside the range of relatively stable nu- 



TABLE II. Dirac-Fock calculation of energy levels for selected 
ions belonging to the ionization sequence of tungsten, com- 
pared with available data. 



Ion 


Ionization Energy 


(10 d cm L ) 


% Deviation 




DF Calc. 


Expt. [35] 




w b + 


-509 


-522 


2.49 


W 11+ 


-1846 


-1868 


1.17 


w 13+ 


-2440 


-2345 


4.05 


w 27 + 


-7075 


-7109 


0.47 


w 37 + 


-13049 


-13080 


0.23 


w 45 + 


-19487 


-19471 


0.08 


w 55+ 


-43101 


-43133 


0.07 


w 73+ 


-652346 


-651338 


0.15 



clei [Z > 100). The most interesting cases are those that 
occur in HCIs with Z ion > 5: crossings c (4/ — 5s), d 
(4/ — bp), and h (5/ — 6p), which are studied in further 
detail in Sections IIV1 IVl and I VII respectively. Once an 
ion is found with orbitals near a particular level cross- 
ing, nearby ions with the same crossing can generally 
be found by increasing the nuclear charge while simulta- 
neously increasing the number of electrons by the same 
amount, provided that the orbital shells involved in the 
crossing are not completely filled. 

a. 3d — 4s: The earliest crossing possible in the 
periodic table occurs in the K isoelectronic sequence 
(N = 19). The ground state configuration is [Kr]4s, but 
the ground state of Sc 2+ (Z = 21) is [Kr]3<£ This cross- 
ing can be seen in the early transition metals, where it is 
well known that the 3d and 4s orbitals are nearly degen- 
erate in neutral and near-neutral ions of these elements. 
All isoelectronic sequences beginning from neutral atoms 
with 19 < N < 28 have this crossing. The N = 29 
isoelectronic sequence starts with Cu, where the ground 
state is 3d 10 4s; this sequence has no crossing since in 
the neutral atom the electron shells already fill in the 
Coulomb-limit order. 

b. 4d — 5s: For the Rb isoelectronic sequence, this 
crossing point occurs near Z = 39, which is Y 2+ . Again 
this level crossing happens in near-neutral systems; it 
is available in isoelectronic sequences with 37 < N < 
46. For N — 47 the ground state already has Coulomb 
degenerate ordering. One ion with this crossing, the two- 
valence-electron ion Zr 2+ , was discussed in [361 ] . 

c. 4/ — 5s: The 5s and 4/ level crossing occurs at a 
higher degree of ionisation than the previous two cross- 
ings. The lightest ions with this crossing occur in the 
N = 47 isoelectronic sequence, which has a single elec- 
tron above closed shells. The ions Nd 13+ , Pm 14+ and 
Sm 15+ have optical transitions between these orbitals; 
they were studied in [27j ■ The heaviest ions with this 
crossing occur when the 4/ and 5s shells are nearly 
filled, i.e. in the isoelectronic sequences of Pm or Nd. 
These were studied in [28| where the ions Ir 16+ and Ir 17+ 
(ground state configurations 4/ 13 5s 2 and 4/ 13 5s, respec- 
tively) were found to have optical transitions from the 
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ground state with the extremely large q values. The to- 
tal number of ions with this crossing is around fifty. This 
level crossing is available in isoelectronic sequences with 
47 < N < 61. We discuss other examples with this cross- 
ing in Sec. HVl 

d. 4/ — bp: The 5pi/2 and 5j>3/2 orbitals are sep- 
arated by a large fine-structure interval, which causes 
this level crossing to occur over a wider range of Z (see 
Fig. [2]). For a single electron above a closed shell, this 
crossing occurs at around Z = 59. Fig. [3] illustrates the 
effect of including MBPT corrections on the position of 
the level crossing. This level crossing is available in iso- 
electronic sequences with 49 < N < 67. The ions W 7+ 
and W 8+ (N = 67 and 66, respectively), which have hole 
transitions between the nearly filled shells, were studied 
in detail in [28|. We discuss other examples in Sec. IVl 

e. 4/ — 6s: This crossing point occurs much ear- 
lier in the ionization sequence than other s — / crossings 
presented here since the difference in principal quantum 
number between the orbitals is An = 2. In the Cs iso- 
electronic sequence, this crossing occurs in Ce 3+ , however 
the 5d orbital also plays a role here and the 4/ — 6s level 
crossing is not seen in the ground state of this sequence. 
This level crossing is available in isoelectronic sequences 
with 55 < N < 69. 

/. 5c? — 6s; Just as in the Ad — 5s case, the 5d and 
6s orbitals cross at low ionization stage; for the Cs iso- 
electronic sequence (N = 55) it occurs in doubly-ionised 
lanthanum [Z = 57). On the other hand s 2 — d 2 tran- 
sitions can have reasonably large g-values even in ions 
with relatively small ion stage, especially where the hole 
transitions are used. Several interesting examples includ- 
ing Hf 2+ , Hg 2+ , and Hg 3+ were studied in [36]. This 
level crossing is available in isoelectronic sequences with 
55 < TV < 78. 

g. 5/ — 6s: This crossing is similar to the 4/ — 5s 
crossing previously discussed, and it was hoped that ions 
which showed this crossing would have very high q- values 
due to the large Z 2 enhancement factor. However, the 6s 
orbital is much more tightly bound than the 5/ orbitals 
and as a result the level crossing occurs at Z = 105 for 
the Au isoelectronic sequence and well beyond 105 for 
the Tl isoelectronic sequence. While this level crossing 
occurs in isoelectronic sequences with 79 < N < 101, it 
is unavailable in any stable nuclei. 

h. 5/ — 6p: The 6p±/2 and 6^3/2 orbitals are very far 
apart in HCIs due to large fine-structure splitting (the 
5/5/2 an d 5/7/2 orbitals are much closer). This causes 
a bifurcation of this level crossing, with 5/ crossing the 
6^3/2 orbitals (in the excited state) near Z = 93 and 
6pi/2 near Z — 98 for the Tl isoelectronic sequence. This 
level crossing is available in isoelectronic sequences with 
81<Z<101. It was originally exploited in [3?} where 
it was shown that optical transitions in Cf 16+ (N = 82 
with two valence electrons) have the largest sensitivity 
to variation of the fine-structure constant seen in any 
atomic system. We discuss more examples in Sec. IVI1 



i. 5/ — 7s: As in the case of the 4/ — 6s crossing, the 
difference in principal quantum number is An = 2. Since 
states with larger n tend to have lower orbital energy, 
this causes the 7s orbital to be comparable in energy to 
5/, thus creating a crossing point early in the ionization 
sequence. Ac 2+ , which is near this level crossing, was 
examined in [3^|. This level crossing is available in iso- 
electronic sequences with 87 < N < 101. 

j. 6d — 7s and 6/ — 7s: The 6c? — 7s level crossing 
occurs in low ionisation stages of isoelectronic sequences 
with N > 87. For example in the Fr isoelectronic se- 
quence, it is seen in Ac 2+ (36|, which has 7p, 6d, and 5/ 
orbitals all within optical range of the 7s ground state. 
This level crossing is available in isoelectronic sequences 
with N > 87. The 6/ — 7s crossing should exist in all se- 
quences with TV > 87 since the 6/ shell is never occupied. 
However the 6/ orbital is at such high energy that the 
crossing occurs in very highly charged ions, with Z > 100. 
Therefore the crossing is not shown in Fig. [T] since it will 
not occur in stable isotopes. 



IV. 4/ - 5s CROSSING 

In this section, we examine the 4/ — 5s crossing in 
greater detail. As mentioned previously, this level cross- 
ing occurs in ions with a relatively high degree of ion- 
ization. Table IIIII presents CI calculations of some ions 
near this crossing with up to three valence electrons. As 
can be seen from the tables, the range of values for the 
charge on the ion Z- lon for which the level crossing occurs 
remains relatively stable. This reinforces the general rule 
of thumb that given an ion near a level crossing, simulta- 
neously increasing or decreasing both the charge Z and 
the number of electrons N by the same amount will result 
in another ion near the same level crossing. 

The 4/ — 5s level crossing is particularly unique in 
that it is the only available level crossing between levels 
of different parity in HCIs. This hints to the possibility 
of optical El transitions in these ions, which could be 
useful for cooling and trapping of HCIs. Two points are 
worth noting, however: firstly, in HCIs the strength of 
El transitions is suppressed compared to near-neutral 
atoms (Sec. IVII B[) ; and secondly with Al = 3 for an 
s — / transition, it will tend to proceed via configura- 
tion mixing which greatly reduces its strength. Exam- 
ples include g3Eu 14+ which has El transitions between 
the ground state, 4/ 2 5s (J p = 3.5 + ), and excited states 
4/5s 2 (J = 2.5") and 4/ 3 (J = 4.5~), with energy inter- 
vals 17478 cm -1 and 28828 cm -1 , respectively. 

Perhaps the most interesting examples presented in 
Table IIIII are the two- valence-electron ion Sm 14+ , which 
was studied in [27J, and the three- valence-electron ion 
Eu 14+ . Both of these ions have ground states with half- 
open 5s shell, which means that both have optical s — f 
and f — s ground state transitions. The two El tran- 
sitions in Eu 14+ mentioned previously are of this type, 
and it means that they will have q- values that are of op- 
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posite sign. On the other hand they may be too broad 
for high-precision clocks. Better reference transitions for- 
elocks are strongly suppressed El transitions, suppressed 
Ml, and E2 transitions. 

It is also possible to have level crossings in hole states, 
where one or two electrons are removed from otherwise 
closed shells and effectively give rise to a similar struc- 
ture as one- or two- valence-electron systems. The specific 
cases of Ir 16+ and Ir 17+ were studied in [28[ for the hole 
case, which leaves all intermediate cases. In general, in- 
termediate ions with more than one electron result in 
large configuration spreading, significantly complicating 
the level structure of the ion. This is not true for hole 
cases, which allow for simpler level structures, yet pro- 
viding the benefit of increased Z and therefore high sen- 
sitivity to a-variation. 

TABLE III: Configuration interaction calculations for the 
level structure of highly charged ions with one, two or three 
valence electrons and 4/ — 5s intervals below 100 000 cm -1 . 
In general, the ellipses . . . used indicate that there are more 
fine-structure states available which we omit for brevity. 



N 


Ion 


Config. 




Energy (cm 1 ) 


47 


TV T .I'J 1 

6oNd ld+ 


5s 


„ __|_ 

0.5+ 









If 


2.5" 


64084 






if 


3.5" 


68480 


47 


6iPm 14+ 


5s 


0.5+ 









if 


2.5" 


8902 






4/ 


3.5" 


14290 


47 


1 r i-l- 

62 Sm 15+ 


4./' 


2.5" 









if 


3.5" 


6485 






5s 


0.5+ 


51314 


48 


eoNd 12+ 


5s 2 


0+ 









4/5s 


2" 


86136 






4/5s 


3" 


87464 






4/5s 


4" 


90435 






4/5s 


3" 


96929 


48 


6iPm 13 + 


5s 2 


0+ 









4/5s 


2" 


32742 






4/5s 


3" 


34261 






4/5s 


4" 


38030 






4/5s 


3" 


44299 






if 


4+ 


98912 


48 


6 2 Sm 14 + 


4/5s 


2" 









4/5s 


3" 


1697 






4/5s 


4" 


6381 






4/ 2 


4+ 


11223 






4/5s 


3" 


12460 






4/ 2 


5+ 


16095 






5* 2 


0+ 


25338 






4/ 2 


4+ 


31069 


48 


6sEu 15 + 


4? 


4+ 









if 


5+ 


5886 






4/5s 


2" 


48780 






4/5s 


3" 


50643 


49 


62Sm 13+ 


4/5s 2 


2.5" 









4/5s 2 


3.5" 


6189 



TABLE III: (continued) 



N 


Ion 


Config. 


J p 


Energy (cm i ) 






if 5s 


3.5+ 


40211 






4/ 5s 


4.0 


A 0/1 £*A 
4Z404 


49 


63EU 14+ 


4/ 2 5s 


3.5+ 









4/ 5s 


A £ + 
4.0 


ZDUl 






4/ 2 5s 


5.5+ 


6663 






4/ 2 5s 


1.5+ 


10711 






4/5s 2 


2.5" 


17478 






4/ 2 5s 


3.5+ 


20722 






4/5s 2 


3.5" 


24854 






A f 3 


A ^ — 
4.0 


ZooZo 


49 


64Gd 15 + 


if 


4.5" 









4/ 


0.0 


4 (Do 






if 


6.5" 


9711 






if 


1.5" 


24137 






if 5s 


3.5+ 


30172 






A f3 


4.0 




49 


B5 Tb 16+ 


if 


4.5" 









if 


5.5" 


5702 






if 


6.5" 


11527 






if 


1.5" 


25637 






4/ 2 5s 


3.5+ 


94034 






4/ 2 5s 


4.5+ 


97331 



V. 4/ - 5p CROSSING 

The 4/ — 5p crossing differs from the 4/ — 5s crossing 
in that the orbitals are of the same parity and the 5p or- 
bital has strong fine-structure splitting. HCIs near this 
crossing can have Ml transitions even without configura- 
tion mixing since the single-electron 5p3/2 ~ 4/ 5 / 2 tran- 
sition is Ml-allowed (although not in the non-relativistic 
limit since Al = 2). Since the ratio of Ml/El transi- 
tion strengths is larger in HCIs relative to near-neutral 
ions (due to the suppression of El transitions), these ions 
can have rich physics to exploit in clocks. Additionally 
E2-allowed transitions are plentiful, and these can have 
linewidths which are more appropriate for reference tran- 
sitions. 

With the possibility of up to 14 electrons in the /-shell 
and 6 electrons in the p-shell, there are many ions that 
have this crossing, from single- valence-electron examples 
like 5gPr 10+ to the nineteen- valence-electron (single hole) 
74 W 7+ . Table [V] presents Dirac-Fock calculations of en- 
ergy levels in the V^" 1 approximation, with a 5p x shell 
included above the closed Cd (N = 48) core. The hp x 
shell (x = N — 1 — 48) is included by weighting the poten- 
tial of the filled 5p 6 shell by the factor x/6. The position 
of the crossing from the Dirac-Fock estimate does not al- 
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TABLE IV. Weighted Dirac-Fock energy intervals calculated 
in the V N_1 potential for highly charged ions near the Af — bp 
level crossing. The Dirac-Fock procedure includes a Cd core 
and a weighted 5p shell: [Kr] bs 2 Ad 10 bp x , with x = N — 49. 



N x 


Ion 


Energy relative to 4/ 5 / 2 


orbital (cm 1 ) 






5pi/ 2 


5^3/2 


4/7/2 


49 


57La B+ 


-114738 


-86257 


1698 




~ q i 

5sCe 9+ 


-70791 


-37066 


2406 




t-. i n_i_ 

59Pr lu+ 


-20261 


19230 


3200 




H.T ill -1- 


36280 


82094 


4085 




6 iPm 12+ 


98419 


151157 


5066 


50 1 


TVT 1 10-1- 


-130828 


-82538 


4013 




i-» 114- 


-77060 


-21675 


4988 




n 124- 

62om 


-17868 


45257 


6065 




eaEu 1 ^ 


46466 


118020 


7249 


51 2 


t-» 1 04- 


-100727 


-47479 


4916 




n 114- 

62Sm 11+ 


-43378 


17463 


5987 




-r-i 1 94- 

63 Eu 12+ 


19191 


88305 


7164 




i 1 ^4- 


86736 


164841 


8456 


52 3 


6 iPm 9+ 


-123164 


-72027 


4849 




1 04- 

62Sm lu+ 


-67686 


-9103 


5915 




-j — . 114- 

63EU 11 + 


-6903 


59797 


7086 




^, ,1 94- 

64 Gd 12+ 


58917 


134449 


8370 


53 4 


6iPm 8+ 


-144342 


-95288 


4787 




i-i Q4- 

62Sm J+ 


-90768 


-34415 


5849 




6 3 Eu 10+ 


-31796 


32520 


7014 




64Gd U + 


32282 


105269 


8290 


54 5 


62 Sm 8 + 


-112596 


-58445 


5788 




6 3 Eu 9 + 


-55462 


6497 


6948 




64Gd 10 + 


6852 


77323 


8218 




65 Tb 11+ 


74101 


153821 


9606 


55 6 


62Sm 7+ 


-133139 


-81160 


5734 




63 Eu 8 + 


-77874 


-18239 


6889 




64 Gd 9+ 


-17346 


50637 


8153 




6 5 Tb 10 + 


48176 


125243 


9533 


ways 


jree with the confij: 


juration interaction calculation. 



but at least provides for a reasonable starting point. 

In Table [V] we present CI calculations for two- and 
three- valence-electron ions near the 4/ — bp crossing. In- 
teresting examples here include the 6oNd 10+ ion which 
has a mixed Afbp ground state from which narrow tran- 
sitions are available to 4/ 2 and bp 2 configurations. These 
would have q values of opposite sign, and so a clock us- 
ing these transition would be a good probe of a- variation. 
On the other hand configuration mixing ensures that the 
three- valence-electron ions 6oNd 9+ and 6iPm 10+ ion have 
good E2 and Ml transitions well within the range of usual 
optical and near-IR lasers. At the heavier end of the 
spectrum of ions which have this crossing are W 7+ and 
W 8+ , with one and two holes in otherwise filled orbitals, 
respectively. These were studied in [28j . 



VI. 5/ - 6p CROSSING 



TABLE V. Configuration interaction calculations for the level 
structure of HCIs with two or three valence electrons and 
4/ — bp intervals below 100 000 cm -1 . The ellipses . . . are used 
to indicate that there are more fine-structure states available 
which we omit for brevity. 



N 
50 



50 



Ion 



50 



51 



3 Pr 8 + 



51 



3 Nd 9 + 



51 



iPm 1 



Config. 


J 


Energy (cm 1 ) 


bp A 








5p 2 


1 


23362 


bp 2 


2 


31033 


4/5p 


3 


92661 


4/5p 


4 


98806 


5p 2 








5p 2 


1 


28273 


5p 2 


2 


34999 


4/5p 


3 


44738 


4/5p 


4 


51669 


4/5p 


5 


86593 


4/5p 


3 





4/5p 


2 


3640 


4/5p 


4 


7701 


bn 2 
op 





9060 


if 


5 


33730 


4/ 2 


6 


36668 


bp 2 


1 


42578 


5p 3 


1.5 





5p 3 


1.5 


26953 


op 


j~ . 


34494 


4/5p 2 


2.5 


47413 


A f ^ 2 
4 J oP 


. 


^0097 


bp 3 


0.5 


51929 


4f5p 2 


3.5 


68470 


4f5p 2 


2.5 





Afbp 2 


3.5 


6429 


bp 3 


1.5 


10613 


Afbp 2 


2.5 


25124 


4/ Op 


a . O 


z / 041 


Afbp 


4.5 


58361 


Afbp 2 


2.5 





Afbp 


4.5 


3937 


Afbp 2 


3.5 


6992 


Afbp 


3.5 


9483 


Afbp 


5.5 


10844 


Afbp 


2.5 


13732 


Afbp 


3.5 


16000 


Afbp 2 


2.5 


31646 



The bf — 6p level crossing is similar in many ways to the 
Af — bp crossing, with some important differences. Since 



this crossing occurs in ions with very high Z, the fine- 
structure splitting of the 6p levels is very large, and near 
the level crossing is usually much larger than the 5/ — 6p 
interval. This provides advantages over the Af — bp cross- 
ing in that there are a larger number of ions available 
where one of these orbitals cross and also that the large 
fine-structure splitting causes a simplification of the level 
structure. In cases where the Qpi/2 and 5/ levels cross, 
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TABLE VI. Weighted Dirac-Fock energy intervals calculated 
in the V N_1 potential for highly charged ions near the 5/ — 6p 
level crossing. The Dirac-Fock procedure includes a Hg core 
(N = 80) and a weighted 6p shell: [Xe] 6s 2 5d 10 4f 14 6p x , with 
x = iV-81. 



N 


Ion 


Energy relative to 5/5/2 


orbital (cm 1 ) 






6P1/2 


6^3/2 


5/7/2 


81 


9oTh a+ 


-182828 


-94604 


4964 


81 


9iPa 10+ 


-167192 


-65878 


6440 


81 


92U 11 + 


-148656 


-33232 


8044 


81 


9 3 Np 12+ 


-127546 


3067 


9775 


81 


94 Pu 13+ 


-104132 


42813 


11636 


81 


9 5 Am 14+ 


-78639 


85844 


13630 


81 


9 6 Cm 15 + 


-51267 


132033 


15760 


81 


97 Bk 16+ 


-22195 


181274 


18030 


81 


9 8 Cf 17+ 


8415 


233481 


20446 


81 


9 9 Es 18 + 


40410 


288583 


23012 


81 


iooFm 19+ 


73643 


346516 


25732 


82 


95Am 13+ 


-258013 


-86547 


14560 


82 


t~i 14+ 

96 Cm 


-237735 


-47179 


16727 


82 


97 Bk 15 + 


-215633 


-4627 


19033 


82 


9 8 Cf 16+ 


-191886 


41006 


21484 


82 


99ES 17+ 


-166663 


89635 


24085 


83 


96Cm 13 + 


-248665 


-63867 


16320 


83 


97 Bk 14 + 


-227341 


-22372 


18614 


83 


9 8 Cf 15 + 


-204330 


22234 


21052 


83 


9 9 Es 16 + 


-179808 


69859 


23638 


84 


96Cm 12 + 


-259040 


-79947 


15920 


84 


97 Bk 13 + 


-238508 


-39525 


18205 


84 


9 8 Cf 14+ 


-216246 


4043 


20630 


84 


99ES 15 + 


-192434 


50657 


23202 


85 


96Cm 11+ 


-262798 


-89202 


15076 


85 


97 Bk 12 + 


-242718 


-49486 


17332 


85 


9 8 Cf 13 + 


-220875 


-6600 


19728 


85 


9 9 Es 14+ 


-197458 


39345 


22269 


85 


iooFm 15 + 


-172644 


88261 


24960 


86 


97 Bk 11+ 


-259178 


-71990 


17413 


86 


9 8 Cf 12 + 


-238450 


-30538 


19816 


86 


9 9 Es 13 + 


-216091 


14017 


22363 


86 


iooFm 14 + 


-192285 


61579 


25059 


87 


97 Bk 10 + 


-268654 


-87268 


17032 


87 


98Cf 11+ 


-248716 


-46899 


19426 


87 


99ES 12 + 


-227101 


-3390 


21960 


87 


iooFm 13+ 


-204002 


43149 


24643 


87 


ioiMd 14+ 


-179597 


92630 


27480 



such as near Cf 17+ , there is an enhancement of sensitiv- 
ity to a- variation [37(. The lower component of the P1/2 
Dirac spinor has an Sj/2 structure and is not small be- 
cause of the high Z. This means that the pi/2 orbital has 
a g-value comparable to an s-wave orbital. 

In Tabic IVII we present weighted Dirac-Fock orbital 
energies for ions near the 5/ — Qp crossing in the V N_1 
approximation. For the single- valence electron case (JV = 
81) two crossings are seen. The first occurs between U 11+ 
and Np 12+ and corresponds to the 5/ — 6p 3 / 2 crossing, 
while the 5/ — Qp\/2 crossing occurs near Cf 17+ . This 
second crossing is only shown in Table IVII for N = 81 



TABLE VII. Configuration interaction estimates for the level 
structure of highly charged ions with two valence electrons 
and 5/ — 6p intervals below 100 000 cm -1 . Ellipses (...) indi- 
cate that there are more fine-structure states that have been 
omitted. All levels have even parity. 



N 



Ion 



82 97 Bk 15+ 



82 98 Cf lfH 



82 



82 



iFm 



18+ 



Config. 


J 


Energy (cm x ) 


Qp A 








5/6p 


3 


89786 


5/6p 


2 


97898 


6p 2 








5/6p 


3 


63664 


5/6p 


2 


72221 


5/6p 


3 


83564 


5/6p 


4 


85846 


6p 2 








5/6p 


3 


36004 


5/6p 


2 


44444 


5/6p 


3 


58033 


5/6p 


4 


59702 


5/ 2 


4 


90788 


6p 2 








5/6p 


3 


7452 


5/6p 


2 


14775 


5/ 2 


4 


28824 


5f6p 


3 


31436 


5/ 2 


4 


36157 


5/ 2 


4 





5/ 2 


2 


4536 


5/6p 


3 


5323 


5/ 2 


5 


20460 


5/ 2 


4 


21371 


5/6p 


2 


22436 


6p 2 





22871 


5/ 


3 


24/95 


5/6p 


3 


35959 


5/ 2 


6 


42232 


5/ 2 


4 





5/ 2 


2 


9828 


5/ 2 


5 


22772 


5/ 2 


4 


27627 


5/ 2 


3 


28876 


5/6p 


3 


37642 


5/ 2 


6 


39897 


5/6p 


3 


67819 


5/6p 


4 


73676 



because it is soon pushed to ions with Z > 100 (although 
clearly it will still occur for two- or three- valence-electron 
ions). 

Configuration interaction calculations for some inter- 
esting HCIs with the 5/ — 6p crossing are shown in Ta- 
ble lVIII (two- valence-electron ions) and Tabic [VIlT] (three- 
valence-electron ions). As with the ions near the 4/ — bp 
level crossing, many Ml and E2 transitions are available 
within the optical range corresponding to single-electron 
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TABLE VIII. Configuration interaction calculations for the 
level structure of HCIs with three valence electrons and 5/ — 
6p intervals below 100 000 cm -1 . Ellipses (. . .) indicate that 
there are more fine-structure states that have been omitted. 
All levels have odd parity. 



N 



Ion 



83 



83 9 6 Cm 13+ 



S3 



r Bk 14+ 



83 



sCf 1 



83 



99ES 1 



83 



jFm 



17+ 



Config. 


J 


Energy (cm 1 ) 


6p d 


1.5 





5/6p 2 


2.5 


70210 


5/6p 2 


3.5 


82945 


5/6p 2 


2.5 





^ f fin 2 

oj op 


o . u 


1 SSI ^ 


5/ 2 6p 


4.5 


83815 


5/ 2 6p 


2.5 


97251 


5/ 2 6p 


3.5 


97765 


5/6p 2 


2.5 





5/6p 2 


3.5 


20858 


5/ 2 6p 


4.5 


58127 


5/ 2 6p 


2.5 


72447 


^ f 2 fin 


O . J 


7^1 8Q 


5/ 2 6p 


1.5 


76551 


5/ 2 6p 


5.5 


79687 


5/6p 2 


2.5 





5/6p 2 


3.5 


22742 


5/ 2 6p 


4.5 


31188 


5/ 2 6p 


2.5 


46699 


^ f 2 fin 
0/ Dp 


O . J 


4-71 ^fi 

*± l lOD 


5/ 2 6p 


1.5 


49751 


0/ bp 


G G 
O.D 


O4oyo 


5/6p 2 


2.5 





o/ op 


4.0 




5/ 2 6p 


3.5 


19106 


5/ 2 6p 


1.5 


22246 


5/ 2 6p 


2.5 


23262 


5/6p 2 


3.5 


26967 


5/ 2 6p 


5.5 


30767 


5? 


5.5 


55606 


5? 


6.5 


64091 


5/ 3 


1.5 


65019 


5/ 2 6p 


4.5 





5/ 3 


4.5 


8162 


5/ 2 6p 


2.5 


11213 


5/ 2 6p 


1.5 


12028 


5f 2 6p 


3.5 


18134 


5/ 3 


5.5 


22763 


5/ 2 6p 


2.5 


28805 


5? 


1.5 


33451 



p — f transitions. Clearly the difficulty with exploiting 
this crossing is that many of the elements with transitions 
near it are not stable and do not occur naturally. In 37] 
we studied Cf 16+ in some detail since it is relatively sta- 
ble (with isotopes that live up to several hundred years) 



and has the &P1/2 — 5/ crossing mentioned previously. 
Using hole transitions is not possible with this crossing 
since there would need to be around 14 or 15 valence 
electrons (corresponding to the crossing of filled 6p 2 / 2 

and 5/ 14 shells, minus one or two electrons), and this 
would require Z > 100, well past the somewhat stable 
elements. 

An interesting example that makes use of the 5/— 6p 3 / 2 
crossing is the three- valence-electron U 9+ . Because of the 
large fine-structure splitting, the first two valence elec- 
trons fill the QPi/2 subshell. The third valence electron is 
in the 6p 3 / 2 subshell (ground state) but may be excited 
to the 5/5/2 and 5/7/2 orbitals. These transitions are 
shown in Table NTU\ The transitions (Ml at 70210 cm" 1 
and E2 at 82945 cm -1 ) will not be particularly sensitive 
to ev-variation, but 235 U is interesting also because it has 
a 76 eV nuclear transition which may soon come within 
XUV laser range [38| . The nuclear transition would have 
high sensitivity to variation of fundamental constants, 
and the electronic transition could then form an 'anchor' 
(relatively insensitive) transition. 



VII. APPROXIMATE SCALING LAWS 

It is useful to have simple estimates of various proper- 
ties of HCIs given existing knowledge of an appropriate 
neutral or near neutral ion. In the following sections we 
derive approximate scaling laws for highly charged ions 
that state how a particular property will change along an 
isoelectronic sequence with increasing Z. Our approach 
is similar to that of [l[ , however we use the formalism of 
effective charges which we believe will be more useful to 
experimentalists. A summary of our results is presented 
in Table KU 



A. Coefficients of Linear Fitting for Effective 
Charge 

Recall the approximate formula for the non-relativistic 
energy of an electron in a screened Coulomb potential, 
V(r) ~ -Z a /r, 



E n 



{Z ion + 1) S 

2i/ 2 



3. 

2n 2 



(3) 



where n is the integer principal quantum number and v 
is the effective principal quantum number which is intro- 
duced to keep agreement with experimentally observed 
energies. For the purposes of this work on highly charged 
ions, it is more convenient introduce an effective charge 
Z a as an alternative to v that represents the (non-integer) 
effective screened charge of the potential that the ex- 
ternal electron 'sees'. In this formulation n is kept as 
the usual integer principal quantum number. Z a scales 
nearly linearly along an isoelectronic sequence as Z (or 
equivalently Z- lon ) increases. 



FIG. 4. Calculated effective charge Z a — ^/\2n' 2 E\ (circles) 
versus ion charge Zi on for a valence 4/ electron above a closed 
shell [Xe] 6s 2 (N — 56) core. The lines represent linear fits 
using the values tabulated in Table IIXI for the appropriate 
regions of Z laTL . 




I i i i i I i i i i I i i i i L_ 

5 10 15 
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Z a can easily be calculated using Dirac-Fock energies 
for any ion. We also present fitting laws for Z a as a func- 
tion of Z; on for several valence orbitals in Table IIXI that 
may be used to quickly obtain Z a . The data were ob- 
tained from one- valence-electron Dirac-Fock calculations, 
and we fit for the linear coefficients A and B according 
to 

Z a = AZ ion + B. (4) 

Values for A and B are presented in Table IIXI For large 
Z ion (5 < Zon < 20, labelled H in Table EX}, our cal- 
culations show that the linear approximation used above 
is in very good agreement with the calculated trend (see 
Fig. [4j. Often experimental data is available for neu- 
tral or near-neutral ions and in order to extrapolate from 
these ions to HCIs requires a reasonable estimate of Z a 
for these ions. Therefore we also present fits across the 
domain 1 < Zi on < 4 (labelled L in Table IIX[) and values 
for the neutral atoms Z- lon = (labelled N in Table IIX|) . 

Our approach is similar to that of Slater [39| for cal- 
culating the effective charge of electrons with shielding 
and the work that followed it. We see that the effective 
charge Z a is always bigger than Z + 1 (this correspond- 
ingly leads to v < n for the other convention). This 
is because electrons spend a non-zero amount of time 
closer to the nucleus, and during that time experience a 
larger ion charge. It is worth noting these two schemes 
are equivalent for the 4/ electron in neutral La, which 
experiences an effective charge very close to Z* + 1 = 1 
because it is much further from the nucleus than the s, p 
and d electrons below it. Therefore the use of Zi + 1 gives 
similar results to Z Q for electrons that are far removed 
from the potential of other electrons. 
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TABLE IX. Table of coefficients A and B for the effective 
charge Z a = AZ\ OTL + B. For the regime column, L means 
the coefficients are more suitable for ions with 1 < Z[ on < 4 
and H means that the coefficients are tailored for high ion 
charge 5 < Zi on < 20, while N is the special case of neutral 
atoms where Z laTL = 0. These values were tabulated using the 
Dirac-Fock energies of singly-occupied electron orbitals above 
closed shells. 



Orbital 


N 


Regime 


A 


B 


4s 


19 


L 


1.233280 


2.4530 






H 


1.060127 


3.3022 






N 




2.1725 


5s 


37 


L 


1.402867 


3.0221 






H 


1.130601 


4.4015 






N 




2.6390 


6.s 


55 


L 


1.542243 


3.4892 






H 


1.198431 


5.2507 






N 




3.0283 


7s 


87 


L 


1.750825 


4.1439 






H 


1.325290 


6.3323 






N 




3.5841 


4p 


31 


L 


1.362347 


2.8994 






H 


1.107311 


4.1724 






N 




2.5049 


hp 


49 


L 


1.525760 


3.5762 






H 


1.179904 


5.3294 






N 




3.0751 


6p 


81 


L 


1.774450 


4.4369 






H 


1.318796 


6.7611 






N 




3.7917 


3d 


21 


L 


1.188580 


4.0655 






H 


1.049841 


4.7758 






N 




3.9083 


id 


39 


L 


1.447306 


2.8779 






H 


1.119397 


4.5226 






N 




2.3191 


5d 


71 


L 


1.724165 


3.4661 






H 


1.227029 


5.9896 






N 




2.6498 


4/ 


57 


L 


1.801169 


2.6428 






H 


1.194715 


5.6147 






N 




1.0082 


5/ 


89 


L 


2.028852 


2.4084 






H 


1.276494 


6.1027 






N 




1.2579 



B. Scaling of EJ and MJ Matrix Elements 

In this section we present analytical estimates for the 
scaling of the EJ and MJ transition matrix elements. 
We use the following formulae to calculate both the non- 
rclativistic electric and magnetic multipole reduced ma- 
trix elements (a relativistic treatment gives the same re- 
sults). In the equations below, we seek to retain only 
the dependence of Z a in the relevant formulae whenever 
possible. The El matrix element is 

{nl\r\n'l') = [ P nl rP n , v dr (5) 
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where the radial wavefunction P n i far away from the core 
electrons is 



P, 



N„ 




-n + Z + 1.2Z + 2, 



2Z„r 



This allows the Z a dependence of the El integral to be 
calculated as 



rP.Pjdr 



(ZaY 



(6) 



The non-relativistic Ml matrix element does not scale 
with charge as it is a function of the angular momenta. 
Therefore, while the El matrix element decreases for a 
decrease in Z a , the Ml matrix element remains constant. 
In comparing highly-charged ions where (large Z a ) and 
near neutral ions (small Z a ), we see that Ml transitions 
can be as strong as El transitions as the latter decreases 
with increasing Z a . A similar treatment was adopted 
in poj . with effective principal quantum number v (la- 
belled n* in their equations) instead of effective charge 
Z a . For higher multipoles one obtains higher powers of 
the Coulomb radius 



(r n ) 



(IB 



(7) 



where as is the Bohr radius, so that the general scaling 
law for EJ and MJ matrix elements is 



m\<lj\Wi) ~ (Z a ) J 



and 



(Ki\\qj M) Wkj) ~ (Z a ) 



1-7 



(8) 



(9) 



In general, E(J+1) matrix elements have the same Z a 
scaling as MJ matrix elements. 



C. Scaling of Polarizability and Blackbody 
Radiation Shift 

The blackbody radiation shift (BBR) for an adiabatic 
system can be calculated using the formula 



SE 



|(831.9 V/mff^-) a (l + ry) (10) 



where ao is the static dipole polarizability and and ?y is 
a small dynamic correction due to the frequency distri- 
bution, which for the purposes of this estimate we will 
disregard. The valence scalar polarizability of an atom 
in a state v can be expressed as a sum over all excited 
intermediate states n allowed by El selection rules. 



a 



3(2 Jv + 1) 



(u||r||n) (n||r||t?) 
E n — E v 



(11) 



TABLE X. Magnetic dipole hyperfirie coefficients A (cal- 
culated in [42]) and their scaling with increasing Z along 
the lithium isoelectronic sequence. Values of Z a were ob- 
ta ined fro m Dirac-Fock calculations using the relation Z a = 
s/\2n 2 E\ where n is the principal quantum number and E is 
the orbital energy in atomic units. The notation | p means to 
use the values in the previous row of the table. 



Isotope 


/ 


91 


A (MHz) [42] 


Z a 


(.A/gi) 


ZZ a 
Zi + 1 


(A/ 9I )\ P 


ZZ a 1 
Zi+T Id 


^Li 


3/2 


2.1709 


399.34 


1.25 






lBe+ 


3/2 


-0.7850 


-625.55 


2.31 


4.35 


4.20 


11 B 2+ 


3/2 


1.7924 


3603.77 


3.33 


2.53 


2.49 


13^3+ 


1/2 


1.4048 


5642.40 


4.35 


2.00 


2.00 




1/2 


-0.5664 


-3973.68 


5.36 


1.74 


1.74 




5/2 


-0.7575 


-8474.13 


6.36 


1.59 


1.59 


19p6+ 


1/2 


5.2578 


88106.93 


7.37 


1.50 


1.50 



We showed in Eq. ([6]) that the reduced matrix element 
(w||r||n) scales simply as 1/Z a . Also, the dependence of 
the non-relativistic energy on Z a is given by Eq. @ to 
be Z\. Therefore all terms in the summation have the 
same dependence on Z a and the total dependence on Z a 
must necessarily be the same. We must then have 



SE 



a 



(12) 



Eq. (|T2j) suggests that in systems with high effective 
charge (large Z a ) such as highly charged ions, the BBR 
shift will be strongly suppressed as compared to neutral 
systems. 



D. Scaling of the Hyperfine structure 

Operators with large negative powers of radius will 
not follow the Coulomb radius scaling, Eq. ([7]), since 
the wavefunction at small distances cannot be described 
by P n i. Instead we must use the approach of Fermi- 
Segre (see, e.g. [4l|) where the normalised squared wave- 
function at the origin ~ Z(Z\ ory + l) 2 /v 3 . Since v = 
n(Zi on + 1)/ Z a , we then come to the following scaling 
law for the hyperfine A coefficient: 



.4 
91 



zzl 



1) 



where we have factored out the nuclear ^-factor gi which 
varies greatly between nuclei. We compare this scaling 
law with experimental data in Table |X] A similar re- 
sult may be derived for the electric quadrupole hyperfine 
constant B. We should also point out that the widths 
of hyperfine transitions will scale as u> 3 ~ A 3 , therefore 
relaxation of hyperfine structure will occur much faster 
in HCIs. 
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TABLE XI. Scaling dependences for HCIs for various sources 
of systematic shifts in optical clocks. 



TABLE XII. Deviations from Madelung filling (usual periodic 
table filling) in neutral atoms. 



2 nd order Stark shift 
Blackbody shift 
2 nd order Zeeman shift 
Electric quadrupole shift 
Fine- structure 
Hyperfine A coefficient 



suppressed 51 

~ Z 2 Zl /(Zion + 1) 
~ ZZ't/iZion + 1) 



The Zeeman shift is sensitive to the specific fine- and 
hyperfine-structure of the transition, but may be suppressed in 
HCIs due to a larger energy denominator. 



VIII. CONCLUSION 

In this paper, we have discussed all level crossings 
available in the periodic table and their characteristics. 
We separately discussed and identified several highly 
charged ions near level crossings and presented estimates 
for the energy intervals in some of these ions. We also 
calculated scaling laws in terms of the effective screened 
charge Z a for transition matrix elements, energy intervals 
(including fine-structure), blackbody radiation and the 
hyperfine shift - these provide a quick and reliable way 
to estimate size of these atomic properties given knowl- 
edge of these properties in a near-neutral ion. In order 
to facilitate these estimates, we have also tabulated em- 
pirical values for Z a for singly occupied electron orbitals 
above closed shells. Our scaling laws predict the BBR 
shifts in HCIs will be strongly suppressed. On the other 
hand, the hyperfine structure is much more important. 

The potential future applications of HCIs as discussed 
in this paper are clear - the strong dependence of transi- 
tions in HCIs on the variation of the fine-structure con- 
stant makes them good candidates for laboratory tests of 
cosmological a-variation. Finally, the existence of level 
crossings leads to the availability of transitions that can 
be excited by optical lasers. This is an experimental ad- 
vantage HCIs have over nuclear clocks which have also 
been proposed to probe the variation of fundamental con- 
stants [2l|, |24j , but require lasers operating in the peta- 
hertz range to excite neutrons or protons. 
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Appendix: Deviations from Madelung Filling 

In some neutral atoms there are deviations from the 
Madelung order of filling; for example, these deviations 



Element 


Actual Filling 


Madelung Filling 


24 Cr 


[Ar] 


3d b 4s 


[Ar] 3d 4 4s 2 


29CU 


[Ar 

L J 


3d 10 4s 


[Ar] 3d 9 4s 2 


4iNb 


\Kt] 

L J 


4d 4 5s 


[Kr] 4d 3 5s 2 


42 Mo 


Kr 

L J 


4d 5 5s 


[Kr] 4d 4 5s 2 


44 Ru 


\Kt] 

L J 


4d 7 5s 


[Kr] 4d 6 5s 2 


45 Rh 


Kr 

L J 


4d 8 5s 


[Kr] 4d 7 5s 2 


46 Pd 


\Kt] 

L J 


4d 9 5s 


[Kr] 4d 8 5s 2 


47Ag 


Kr 

L J 


4d 10 5s 


[Kr] 4d 9 5s 2 


57La 


[Xe] 

L J 


5d6s 2 


[Xe] 4/6s 2 


58 Ce 


[Xc] 


4/5d6s 2 


[Xe] 4/ 2 6s 2 


64 Gd 


[Xe] 


4/ 7 5d6s 2 


[Xe] 4/ 8 6s 2 


78 Pt 


[Xc] 


4/ 14 5d 9 6s 


[Xe] 4/ 14 5d 8 6s 2 


79AU 


[Xe] 


4/ 14 5d 10 6s 


[Xe] 4/ 14 5d 9 6s 2 


sgAc 


[Rn] 


6d7s 2 


[Rn] 5/7s 2 


9 Th 


[Rn] 


6d 2 7s 2 


[Rn] 5/ 2 7s 2 


9i Pa 


[Rn] 


5/ 2 6d7s 2 


[Rn] 5/ 3 7s 2 


92U 


[Rn] 


5/ 3 6d7s 2 


[Rn] 5/ 4 7s 2 


93 Np 


[Rn] 


5/ 4 6d7s 2 


[Rn] 5/ 5 7s 2 


96 Cm 


[Rn] 


5/ 7 6d7s 2 


[Rn] 5/ 8 7s 2 



are commonly observed in the transition elements. Sim- 
ilar to our treatment of the level crossings due to the 
Coulomb degeneracy, we now examine all available devi- 
ations from the Madelung filling order and characterise 
them. Recall that for level crossings that occur due to 
Coulomb degeneracy, we had crossings of type s — f, s — d 
and p — f only. In contrast, we find only find s — d and 
d — f type crossings in isoelectronic sequences starting 
from atoms that have deviations (see Table IXIII for an 
exhaustive list). 

A specific example is lanthanum, with a ground state 
of [Xe] 5d6s 2 . Here the 5d orbital is filled before the 4/ 
orbital, while from the Madelung rule we would expect 
the 4/ to be filled first. Because 5d has higher n than 4/, 
further along the isoelectronic sequence Coulomb degen- 
eracy will cause the 4/ orbital to be lower in energy than 
the 5d orbital. In Cc + the ground state is [Xe]4/5d 2 . 
Pr 2+ has a ground state of [Xe]4/ 3 which shows that 
all crossings have occurred by Z- lon = 2. We find that 
all crossings caused by the deviation from the Madelung 
rule occur at low ion charge. 

a. 3d — 4s: This crossing occurs due to the addi- 
tional stability offered by half-filled and filled 3d orbitals 
in Cr and Cu respectively. The half or complete filling of 
the 3d orbitals is preferred to a filled 4s orbital, as a re- 
sult one of the 4s electrons in these atoms fill a 3d orbital 
instead. For Cr isoelectronic sequence, the remaining 4s 
electron eventually fills a 3d state instead. 

b. 4d — 5s: In the ground states of Nb, Mo, Ru, Rh, 
Pd and Ag, the 4d shell fills before the 5s shell is closed. 
This is consistent with the calculations done on the Rb 
isoelectronic sequence, which reveals that the 4d and 5s 
orbitals cross at Z = 39. 
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c. Af — 5d: Due to angular momentum and parity 
considerations, there exists optical El transitions in neu- 
tral La, Ce and Gd, as well as for near neutral ions in 
the vicinity of these atoms. The 4/ and 5d orbitals must 
necessarily cross due to Coulomb degeneracy. Our calcu- 
lations show that this crossing occurs at Z = 58 for the 
La isoelectronic sequence. 

d. 5d — 6s: The ground states of Pt and Au show a 
deviation from Madelung filling. According to our calcu- 
lations, the 5d orbital lies below the 6s orbital for Z = 56, 
the next ion in the isoelectronic sequence. The crossing 



occurs near the beginning of the isoelectronic sequence 
because the orbitals are very close in energy to begin 
with. 

e. 5/ - 6d: The ground states of Ac, Th, Pa, U, Np 
and Cm have a single electron in the 6c? 3 / 2 orbital. In the 
example of neutral Th the 6d orbitals lie below the 5/ 
orbitals, but in singly ionized Pa + , the level crossing has 
occured [43j and the 5/ orbitals lie around ~ 5000 cm -1 
below the 6d orbitals. This crossing is also present in the 
Th 3+ ion that has several potential atomic clock transi- 
tions with enhanced sensitivity to a- variation [20| . 
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